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Let T(m,n) denote the number of ways to tile an m-by-n rectangle with dominos. For any fixed m, the 
numbers T(m,n) satisfy a linear recurrence relation, and so may be extrapolated to negative values of n; 
these extrapolated values satisfy the relation T(m, —2 — n) = e m „T(m, n), where E mn = — 1 if m = 2 (mod 
4) and n is odd and where £„, „ = +1 otherwise. This is equivalent to a fact demonstrated by Stanley 
using algebraic methods. Here I give a proof that provides, among other things, a uniform combinatorial 
interpretation of T(m,n) that applies regardless of the sign of 



1. INTRODUCTION 

It has undoubtedly been observed many times that the Fi- 
bonacci sequence, when run backwards as well as forwards 
from its initial conditions, yields a doubly-infinite sequence 
that is symmetrical modulo some minus-signs: 

...,-8,5,-3,2,-1,1,0,1,1,2,3,5,8,... 

We will see that this is just one instance of a more general 
phenomenon involving domino tilings of rectangles. For pos- 
itive integers m and n, let T(m,n) denote the number of ways 
to cover an m-by-n rectangle with l-by-2 rectangles ("domi- 
nos") with pairwise disjoint interiors. When m — 2, the val- 
ues T(m,n) form the Fibonacci sequence; more generally, for 
each fixed m, the values of T(m,n) form a sequence satisfying 
a higher-order linear recurrence relation (see section Q for an 
explanation of why such a recurrence relation must exist). 

These recurrences allow one to extrapolate T(m,n) to neg- 
ative values of n in a natural way, for each fixed positive inte- 
ger m. There are in fact infinitely many recurrences one could 
employ; but, as section [I] explains, the extrapolated values of 
T(m,n) do not depend on which of the recurrence relations is 
used. 

One finds that as long as m is not congruent to 2 (mod 4), the 
sequence of values T(m,n) (with n varying over the positive 
and negative integers) exhibits palindromicity: 

T(m, —2 —n) = T(m,n) for all n. 

When m is congruent to 2 (mod 4), a modified version of this 
holds: 

T(m,-2-n) = (-l) n T(m,n) for all n. 

These properties of the sequence of values T(m,n) can also 
be expressed in terms of the generating function F m (x) — 



L~ = o T(m,n)x", with T(m,Q) taken to be 1. This generat- 
ing function is a rational function of x, as a consequence of 
the recurrence relation governing the values T(m,n). Denote 
this rational function by P m (x)/Q m (x)- Note that the polyno- 
mial Q m (x) encodes the recurrence relation that governs the 
sequence of coefficients T(m,n), at least eventually; if the de- 
gree of P m (x) is less than the degree of Q m (x), then the recur- 
rence relation applies to all the coefficients of the generating 
function, not just ones with n sufficiently large. 

Stanley showed [^] (see parts (d) and (e) of his main The- 
orem) that for m > 1, the degree of Q m is 2 more than the 
degree of P m , Q m (\/x) — Q m (x) when m = 2 (mod 4) and 
Q m (1 /x) = —Q m (x) when m^2 (mod 4). These reciprocity 
relations, first observed empirically by Klarner and Pollack 
[j^], can be shown to be equivalent to the claims made earlier 
about the extrapolated values of T(m,n). 

Stanley's proof makes use of algebraic methods; in particu- 
lar, it relies upon formulas proved by Kasteleyn [^] and Tem- 
perley and Fisher [^[j that express T(m,n) in terms of complex 
roots of unity (or, if one prefers, trigonometric functions). 

I will present a combinatorial way of thinking about the 
two-sided infinite sequence (r(m,n))" = _ 00 . All features of the 
reciprocity relation — including the factor e mn — will have 
direct combinatorial interpretations. Moreover, the method 
used here can also be used to demonstrate similar reciprocity 
phenomena for three- and higher-dimensional analogues of 
domino tiling. For instance, if T(k,m,n) denotes the number 
of ways to tile a k-by-m-by-n box with 1 -by- l-by-2 blocks (in 
any orientation), then T(k,m,— 2 — n) = ±T(k,m,n), where 
T(k,m, —2 — n) is obtained by taking the recurrence relation 
for T(k,m,n) (with k,m fixed and « varying) and running it 
backward. Such higher-dimensional tiling problems are not 
amenable to Stanley's method, since the formulas of Kaste- 
leyn and Temperley and Fisher do not extend beyond the two- 
dimensional setting. 

2. SIGNED MATCHINGS 
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Instead of using domino-tilings, we will use an equivalent 
mathematical model: perfect matchings of grid-graphs, also 



2 



known as dimer covers or dimer configurations 1 . A perfect 
matching of a graph £f (hereafter called just a matching for 
short) is a subset of the edge-set of §f such that each vertex of 
S? belongs to one and only one edge of the subset. The domino 
tilings of an m-by-n rectangle are in obvious bijection with 
the perfect matchings of a graph with mn vertices arranged in 
m rows of n vertices each, with edges joining horizontal and 
vertical neighbors. See Figure 1. 



Fig. 1: The graph ^(2,5). 

Given that some of the extrapolated numbers T(m,n) are 
negative, it is natural to try to provide them with a combi- 
natorial interpretation by way of a signed version of graph- 
matching. A signed graph is a graph whose edges are labeled 
+ 1 and —1, and a matching counts as either positive or nega- 
tive according to whether the number of — 1 edges is even or 
odd. (This is a special case of weighted matching, in which 
each edge is assigned a weight, and the weight of a matching 
is the product of the weights of its edges.) When we count 
the matchings of a signed graph, we count them according to 
sign. For example, the Fibonacci number T(2, —5) = —3 will 
arise as the "number" of signed matchings of the signed graph 
shown in Figure 2. 
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Fig. 2: The signed graph Sf (2, —5). 

Here the horizontal edges are all +1 edges and the vertical 
edges are — 1 edges (where the presence of a negative edge is 
indicated by use of dashed lines). Each of the three matchings 
of the graph has an odd number of vertical edges and so counts 
as a negative matching. Hereafter, when I refer to the number 
of matchings of a signed graph, I always mean the number of 
positive matchings minus the number of negative matchings. 

By a "signed graph of width m", I will mean a subgraph 
of the m-by-n grid graph for some n, with some edges having 
sign 1 and others having sign — 1 . We will be pay special at- 
tention to some specific graphs @(m,n), defined form > 1 and 
n an arbitrary integer. When n is positive, @(m,n) is just the 
m-by-n grid graph described above, with all edges having sign 
1. When n is less than or equal to zero, @(m,n) is a modified 
version of the graph Sf (m,2 — n) in which the rightmost and 
leftmost flanks of vertical edges have been removed, and all 
remaining vertical edges are decreed to have sign — 1 . Thus 
Figure 1 shows the graph Sf (2,5) while Figure 2 shows the 
graph Sf(2, -5). 



The difference between the titles of and [H], compared with the dif- 
ference in the names of the journals in which the two articles appear, clarifies 
the difference between dimers and dominos: Studying dimers configurations 
is applied mathematics. Studying tilings by dominos is not. 



It will emerge that the (signed) number of matchings of 
Sf(m,«) is T(m,n). Note that for all m, Sf(m,0) and Sf (/«, —2) 
each have a unique matching, whose sign is +1. This verifies 
the reciprocity theorem in the case n = (and the case n = — 1 
is trivial). To prove the claim for other values of n (it suffices 
to consider n > 0), we must do more work. 



3. ADJUNCTION OF GRAPHS 

The two definitions of Sf (m,n) given above (one for n > 0, 
the other for n < 0) are not just two distinct combinatorial 
notions, patched together at 0; they fit "seamlessly". To make 
this clearer, we will situate the graphs Sf (m,n) in a broader 
context, in which they function as building blocks. Suppose 
Sfj and Sf, are signed graphs of width m. Define the "adjoined 
graph" Sf[Sf 2 as the signed graph we get when we place S^ 
to the left of Sf 2 and join the rightmost vertices of Sfj to the 
respective leftmost vertices of Sf, using n edges of type +1. 
Figure 3 shows the graph Sf (2,3)Sf (2, —3). 
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Fig. 3: The signed graph Sf(2,3)Sf(2, -3). 

Clearly adjunction is an associative operation, so a product of 
three or more signed graphs of width m is well-defined. 

Let M(Sf ) denote the (signed) number of matchings of a 
signed graph Sf . I can now state the precise sense in which the 
two parts of the definition of Sf (m,n) fit together: 

M(^(m,n^(m,n 2 ) ■ ■■&(m,n k )) 

= M(&(m,n 1 +« 2 H \-n k )) 

for all n v . . . ,n k , regardless of sign. 

To indicate why the preceding equation is true, we will 
prove the case k = 2, which contains all the essential ideas. 

When n l and n 2 are positive, we have Sf '(m,n,)S^ '(m,n 2 ) 
literally equal to the graph Sf (ra,Hj + n 2 ) by virtue of the def- 
inition of adjunction, so this case requires no proof. 

When «j and n 2 are both non-positive, we get situations like 
the one shown in Figure 4, depicting Sf (2, — 2)Sf (2, —3). 
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Fig. 4: The signed graph Sf (2, -2)^(2, -3). 

It is a well-known lemma in the theory of matchings that the 
number of matchings of a graph is unaffected if a chain 
of vertices u,v,w (where v is connected to u and w but to no 
other vertices in the graph) is shrunk down to a single vertex 
(i.e., v is removed and u and w are replaced by a single ver- 
tex adjacent to all the vertices other than v to which u and w 
were originally connected). The same observation applies to 
signed graphs as well if the edges uv and vw are +1 edges. 
Applying this lemma to ^(2, —2)^(2, —3) twice (once in the 
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top row and once in the bottom row), we obtain a copy of the 
graph ^(2, -5). More generally, M(&(m, —a)^(m, —b)) = 
MC$(m,-{a + b))) whenevera,/? > 0. 

When «j and n 2 are of opposite type (one positive, the other 
non-positive), we get situations like the one shown in Figure 3, 
depicting Sf(2,3)£f(2, —3). If we apply the shrinking lemma 
in this case, once in each row, we get the signed graph shown 
in Figure 5. 
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Fig. 5: Sf (2,3)^(2,-3), after shrinking. 

Here the edge marked H — is actually two edges, one a +1 
edge and the other a — 1 edge. Any signed matching that uses 
the + 1 edge can be paired with a signed matching that uses 
the — 1 edge; these matchings have opposite sign, and so to- 
gether contribute to the signed number of matchings. Hence 
the number of matchings of the graph is unaffected if we re- 
move both of these edges, obtaining the signed graph shown 
in Figure 6. 
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Fig. 6: Sf(2,3)Sf(2, -3), reduced to Sf(2,2)Sf(2, -2). 

But this is just §?(2,2)§f(2,— 2). More generally, applying 
the "shrinking lemma" to ^(m,a)^(m, — b) yields &(m,a — 
l)&(m,—(b— 1)) as long as a,b are both positive. This 
sets the stage for a proof by induction, and all that remains 
is to verify the base case where a = or b = 0. This is 
an easy verification that I leave to the reader. The case of 
&(m,—a)^{m,b) is identical. 



4. RECURRENCE RELATIONS 

We start by showing that for every width-m signed graph 
iff, the sequence M(J$f&(m,n)) (with n = 1,2,3,...) satis- 
fies a linear recurrence relation of degree at most 2"' that de- 
pends only on m, not on Jff, 

Let be a width-m signed graph. For any sequence 
Cj,...,c m of m bits, and for n > 1, let S(n;c) be the sum 
of the signs of the matchings of the signed graph obtained 
from (m,ri) by deleting the subset of the rightmost m ver- 
tices specified by the bit-pattern c; specifically, we remove the 
rightmost vertex in row i (and all edges incident with it) for 
precisely those i with c i = 0. Let S(n) be the row-vector of 
length 2 m whose components are the numbers S(n;c). Then 
(cf. [(§ and there is a 2 m -by-2 m matrix M such that 
S(n + 1) = S(n)M. More specifically, the c,c' entry of M 
gives the number of ways in which each matching counted by 
S(n;c) can be extended to a matching counted by S(n + l;c'), 
adding only horizontal edges between the nth and n + 1st 
columns and vertical edges in the n + 1st column. Hence we 
have S(n) = S(l)Af* -1 . By the Cayley-Hamilton Theorem, 



the sequence of powers of M satisfies the linear recurrence 
associated with the characteristic polynomial of M, and it 
follows that the sequence of vectors S(l),S(l)M,S(l)M 2 , . . . 
also satisfies this recurrence. In particular, the entry associ- 
ated with c = (1,1,... ,1) satisfies this recurrence. 
For instance, with m = 2, M is the 4-by-4 matrix 
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10 

10 
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where rows and columns are indexed in the order 
(0,0), (0, 1), (1,0), (1,1). The characteristic polynomial of 
this matrix is x 4 — x 3 — 2x 2 +x+ 1, so the one-sided sequence 
of numbers a n = M(Jf?'&(m,n)) (with n going from 1 to °°) 



-2a n+2 + a n+l +a„ 



must satisfy the recurrence a n+4 — a n+3 
0, regardless of the nature of the width-m signed graph Jtif. 

I now claim that the two-sided sequence of numbers a„ 
(with n going from — °° to °°) satisfies this same recur- 
rence relation. If we take = @(m,—N) with large, 
and use the fact that M(@(m,—N)@(m,n)) = M(&(m,n — 
N)) = T(m,n —N), we see that if we take the doubly-infinite 
sequence . ..,T(m, — l),T(m,0),T(m, 1),... and start it at 
T(m, —N), we get a singly-infinite sequence that (as we have 
shown) must satisfy the Cayley-Hamilton recurrence. Since 
was arbitrary, the two-sided infinite sequence of a n 's satisfies 
the Cayley-Hamilton recurrence as well. 

Having constructed one recurrence relation 
satisfied by the doubly-infinite sequence v = 
(...,T(m,—l),T(m,Q),T(m,l),...), or equivalently, one 
linear operator A that annihilates v, we must now consider 
others. Could there be a linear operator B that annihilates the 
singly-infinite sequence T(m, l),T(m,2),T(m,3),... but not 
the doubly-infinite sequence v? If so, then B would yield a 
different extrapolation of T(m,n) to negative values of n than 
A would. 

To show that no such operator B exists, note that BA annihi- 
lates v. Since BA = AB, B must send v into something annihi- 
lated by A. If B does not itself annihilate v, then B sends v into 
some other sequence that vanishes for all sufficiently large in- 
dices but does not vanish for all indices. It is easy to show that 
such a sequence cannot be annihilated by A or indeed by any 
linear operator with constant coefficients. 

The above arguments show that the numbers T(m,n), de- 
fined for non-positive values of n as above, are the unique 
way to extrapolate so that a linear recurrence is satisfied. 



5. RECIPROCITY 

With all the ingredients in place, we can now give a very 
simple explanation of the reciprocity relation for domino 
tilings of rectangles of fixed width. For concreteness, we start 
with the special case T(2, —5). Referring to Figure 2, we see 
that every matching of Sf (2, —5) must contain the two leftmost 
horizontal edges and the two rightmost horizontal edges. If we 
remove these edges from the graph (and the vertices incident 
with those edges, and the edges incident with those vertices), 
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we get a copy of (2, 3) in which the vertical edges have sign 
1 . Since every matching of J#(2, 3) has an odd number of ver- 
tical edges, we see that T(2, —5) = —T(2,3). 

More generally, when n is positive, every matching of 
( S(m 1 —2 — n) must contain the m leftmost horizontal edges 
and the m rightmost horizontal edges. Removing these edges 
(and concomitant vertices and edges) leaves a copy of &(m, n) 
in which the vertical edges have sign 1 . It is well-known (see 
e.g. [[$]]) that every domino-tiling of a rectangle (or indeed 
more general regions) can be obtained from every other by 
means of moves in which a 2-by-2 block of horizontal domi- 
nos is rotated to give a 2-by-2 block of vertical dominos, or 
vice versa. Since moves of this kind do not affect the parity 
of the number of vertical dominos (or, in our terms, the par- 
ity of the number of vertical edges in a matching), we know 
that that all the signed matchings of our modified version of 

(m, n) carry the same sign; it remains to determine what this 
sign is. 

If one is willing to appeal to the invariance-of-parity result 
mentioned in the preceding paragraph, it is simple to eval- 
uate the common sign of all the matchings by considering 
the all-horizontal matching or the all-vertical matching. But 
in the interests of making the article self-contained, I give a 
direct argument (of a fairly standard kind). Say we have a 
matching of @(m,n) that involves k i vertical edges joining 
row i and row i + 1, for i ranging from 1 to m — 1 . When 
these edges and their vertices are removed from the graph, 
along with all other edges joining row i and row z' + l, our 
graph splits into two subgraphs, each of which must have 
an even number of vertices (since we are assuming that we 
have a matching). Hence ni — k i and n(m — i) — k i must be 
even. Using the congruence fc. = ni (mod 2), we find that 
k x H \-k i, the number of vertical edges, must be con- 
gruent to n(l +2 H h (m— 1)) = n(m — l)m/2 (mod 2). If 

m is congruent to or 1 (mod 4), (m — l)m/2 is even, so the 
number of vertical edges is even. If m is congruent to 3 (mod 
4), then (in order for a matching to exist) n must be even, so 
that n(m — l)m/2 is even. However, if m is congruent to 2 
(mod 4), then (m — l)m/2 is odd, so that n(m — \ )m/2 is even 
or odd according to the parity of n. 

Hence, if we define e m n to be — 1 when m is congruent to 2 
(mod 4) and n is odd, and +1 otherwise, we have 

T(m, —2 — n) = e mtn T(m,n) 

as claimed. 



6. MOTIVATION 

The definition of Sf(m,n) for n < arose not by ad hoc in- 
sight but by a particular vision of combinatorics as a special- 
ized form of algebra. Under this vision, a collection of com- 
binatorial objects ought to be represented by a multivariate 
polynomial whose coefficients are all equal to 1 , where the in- 
dividual terms represent the combinatorial objects themselves 
in some fashion. 

In the case of matchings of ^(2,n) with n > 0, we assign a 
formal variable (or "weight") to each edge, define the weight 



of a matching as the product of the weights of its constituent 
edges, and define the polynomial M(5f(m,n)) as the sum of 
the weights of all the matchings. This is a polynomial in 
which every coefficient equals 1, and in which the constituent 
monomials encode the respective matchings of the graph. If 
we think of Sf(2,n) as being embedded in Sf(2,n + 1) as the 
induced subgraph on the leftmost n columns, we find that the 
edge-variables occurring in M(Sf(2,n)) form a subset of the 
edge-variables occurring in M(Sf (2,n + 1)), and that indeed, 
there is a recurrence expressing M(5f(2,«+ 1)) in terms of 
M(Sf(2,«)) andM(Sf(2,n-l)): 

M(^(2,«+l))=y„ +1 M(^(2,«)) + w„x„M(^(2,«-l)), 

where y n+l is the weight of the rightmost vertical edge of 
&(2,n+ 1) and w n and x n are the weights of the rightmost 
horizontal edges of (2 , n + 1 ) . 

We can run this recurrence backward, obtaining in succes- 
sion some rational functions of the variables w„,x„,y„ with n 
negative. These rational functions are in fact Laurent poly- 
nomials; moreover, all the coefficients of these Laurent poly- 
nomials are ±1. Just as the monomials of M(Sf(2,n)) with 
n > encode the matchings of £f(2,«), one might assume 
that the monomials in the extrapolated Laurent polynomials 
"M(Sf(2,n))" with n < encode combinatorial objects that 
will yield a combinatorial interpretation of the extrapolated 
values T(2,n) with n < 0. The monomials lead one to match- 
ings of the graphs ^(2.n) with n < 0, and once one has 
this signed-matching model, the idea is readily extended to 
T(m,n) for all m. 

It is worth remarking that the reciprocity theorem satisfied 
by the Laurent polynomials M(^(m, n)) is a bit more compli- 
cated than the reciprocity theorem satisfied by the numbers 
T(m,n); in addition to the sign-factor £„,,„, there is also a 
monomial factor involving the variables. A comprehensive 
reciprocity theorem for domino tilings of rectangles would 
take this into account. I intend to write a longer article that 
treats these issues, and also applies them to other models, such 
as matchings of honeycomb hexagons (see [[]]]). I expect that 
algebraic recurrences for sequences of Laurent polynomials 
will also shed light on the issues raised in Problem 32 (the last 
problem) in [jjj. 
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